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ABSTRACT 



The optimal order quantity for the case of uncertainty 
in the newsboy problem depends on the maximum demand. In 
previous work it was assumed that the demand distribution 
was known. In practice, this is often not the case. This 
study suggests some procedures which can be used to estimate 
the demand distribution even if data on unsatisfied demands 
are not available. 
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INTRODUCTION 



I . 



The perishable product is beginning to receive attention 
in the inventory problem. Often, it has a fixed lifetime 
after which it becomes useless to satisfy demands. News- 
papers, medical drugs, human blood, photographic films and 
some foodstuffs are common examples. 

An inventory problem for perishable products is often 
treated as a single-period inventory problem, sometimes 
called a newsboy problem. The newsboy problem was introduced 
by Morse and Kimball in 1950 and the development of it has 
been presented in several journals and publications. A 
brief history of it was given by Masuda [Ref. 1]. 

A well-known characterization of the newsboy problem is : 

The newsboy must decide how many papers to purchase 
for resale during the day. If he buys too many papers — 
that is, more papers than people will purchase from him — 
he will have papers left at the end of the day and will 
incur a loss on each, since old newspapers have little 
salvage value. If, on the other hand, he purchases too 
few papers for resale, he will sell out early and his 
later customers will not be able to buy from him; he thus 
incurs a cost (lost profit) for each paper he could 
have sold, but didn't. The newsboy's problem is to de- 
cide how many papers to have on hand at the beginning of 
the day so as to minimize the losses (due to surplus or 
shortage) at the end of the day. [Ref. 2] 

Finding good solutions to the newsboy problem depends 
upon the amount of information about the magnitude of demand 



X. 



In recent research dealing with newsboy ordering deci- 
sions where little information existed about demand, it was 
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assumed that at the end of each period the demand for that 
period would be known [Refs. 1,3,4]. In order for this to 
happen, one would need unlimited supply or some way to 
measure the demand from unsatisfied customers. In practice, 
this is often not the case. For example, when the newsstand 
is empty, the later customers will not stop at that news- 
stand. Accordingly, we do not know how many more customers 
would like to buy that newspaper. Thus, we may have frequency 
data on demand only up to the value of supply. 

This thesis is intended to provide methods permitting 
solutions to the newsboy problem when we only have demand 
frequencies up to the supply value. 

In Chapter II, two procedures will be proposed to esti- 
mate the maximum demand. The first procedure uses the Gompertz 
Curve since the nature of the demand cumulative distribution 
function (c-d-f.) can be viewed as a growth curve with the 
market saturate as an asymptote. The second procedure sug- 
gested is the four-parameter Schmeiser-Deutsch Distribution 
since it can take on a wide variety of shapes so as to match 
many distributions. Also, the Schmeiser-Deutsch Distribution 
was found to be a suitable procedure for the case of limited 
supply. 

Using the estimation of maximum demand together with the 
unit cost of surplus and the unit cost of outage, the estima- 
tion of optimal order quantity is discussed in Chapter III 
for the "under uncertainty" and "under risk" cases. 



10 



Conclusions and applications are given in Chapter IV. 

A review of the Schmeiser-Deutsch Distribution is presented 
in Appendix A. 
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II. ESTIMATION OF MAXIMUM DEMAND 



Care should be taken to estimate the maximum demand since 
the optimal ordering quantity for uncertainty depends on it. 
In this chapter, two alternatives are suggested to estimate 
the maximum demand when we have unlimited supply. Then, we 
will show how one of these methods may be used to estimate 
maximum demand when there is no data on unsatisfied demand 
(the case of limited supply) . 

The procedures that appear suitable and realistic are 
(1) using the Gompertz Curve and (2) using the Schmeiser- 
Deutsch Distribution. We will discuss the Gompertz Curve 
first. 

A. USING THE GOMPERTZ CURVE TO ESTIMATE THE MAXIMUM DEMAND 
In the terminology of economics, the demand c.d.f. can 
be considered as a "growth curve" since it has a "market 
saturate" as an upper asymptote. The Gompertz Curve is the 
most generally used growth curve. Its formula 

h X 

Y = ka D , (1) 

where X is the time series and Y is the sales amount, de- 
pends on the values of the parameters k, a and b (or their 
logarithms) , it may take on any one of a variety of shapes 
[Ref. 5]. 



12 



Y 




Figure 1. The Gompertz Curve with 0<a<l / 0<b<l 



We are interested in using the Gompertz Curve with 0 < a < 1 
and 0 < b < 1 to approximate the c.d.f. of demand. This form 
of the Gompertz Curve is known to be asymptotic to k. When 
we fit the Gompertz Curve to the observed cumulative fre- 
quency of the data, k should theoretically be 1.0. In prac- 
tice, however, it may depart slightly from 1.0. To compensate 
for this, we use the following form to approximate the c.d.f.: 



Y 

k 



b 

a 



X 



The following description is taken from [Ref. 5] : 

The parameters k, a and b are obtained as follows: 

1) Be certain that the number of observations is divisible 
by 3, i.e., there are 3N periods of base data. 

2) Take the common logarithms of each Y. 

3) Add the first N log Y's to obtain E^ log Y; then the 
second N data points to obtain E 2 log and the last N 
points to obtain E ^ log Y. 

4) Set the smallest observation to be unity. 
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5) Substitute in the following formulas: 



1 3 log y - I 2 log Y 
I 2 log Y - Z log Y ' 



log a 



(Z 2 log Y - ^log Y) • 



(b -1) 

FI 2 
(b -1) ^ 



(3) 



and 



log k = ^ • (It log Y - ( b - — ) • log a). (4) 

N 1 b-1 



With values of k, a and b, one can compute Y for any 
value of X using the formula of the Gompertz Curve. Also, 
the formula can be derived from the Gompertz Curve to 
generate X from the Gompertz Curve as : 



X 



log ( 



log Y - log k . 

log a 

log b 



(5) 



In this chapter, we wish to estimate the maximum demand. 
If the demand distribution has an infinite upper tail, we 
can not hope to estimate the 10 0^ percentile of demand. 

x. I- 

We will be content, therefore, with the 99 percentile of 
demand. This is found by setting (Y/k) = .99 or Y = ,99k 
in equation (5), i.e.. 



X will be 
max 



log ( 



-.0044 
log a 
log b 



( 6 ) 
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which may be considered to be the estimator of the maximum 
demand if the smallest observed X is unity. (If not, the 
maximum demand should be adjusted to be X plus the smallest 
X minus one . ) 

As an example assume that the demands for newspapers 
have been 10, 6, 9, 7, 5, 13, 11, 7, 8, 8. Using the Gompertz 
Curve to estimate the maximum demand, we can preceded as 
follows: 

Let X be the demand for newspapers , and Y be the cumula- 
tive relative frequency of X. 

TABLE I 

DEMAND AND ITS CUMULATIVE RELATIVE FREQUENCY 
X Y 



5 


• — 1 
• 

o 


6 


0.2 


7 


0.4 


8 


0.6 


9 


0.7 


10 


0.8 


11 


0.9 


12 


0.9 


13 


1 . 0 



Since Y is a growth curve, it can be approximated by the 

b X 

Gompertz Curve, Y = ka D . To determine the appropriate 
values of k, a and b using the method described above: 

1) From these data, there are 9 distinct values between 
the smallest observation (5) and the largest observation (13), 
and 9 is divisible by 3. (Here, N = 9/3 = 3.) 
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2) In the third column of Table II, Y's are converted to 



log Y's. 

3) Values of I^log Y, Z 2 log Y and Y are s ^ own i- n 

the fourth column of Table II. 

4) Since the smallest observed X is not 1.0, the maximum 
demand should be adjusted by adding the difference between the 
smallest observed X and 1.0. This is because the Gompertz 
Curve assumes the smallest observation is unity. 

5) From equation (2) ; 



b 



N 



( (-. 0915) -(-.47 37) )/( (-. 47 37) -(-2.0969)) 



0.2355 . 



Therefore 



b 



0.6175 . 



From equation (3); 



log a = (-. 4737-(-2. 0969) ) (.6175-l)/(. 2355-1) 



-1.0622 . 



Therefore 



a 



0.0867 . 



From equation (4) ; 



log k = |(-2. 0969-(. 2355-1) (-1.0622)/(. 675-1) ) 



0.0087. 
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Therefore , 



k = 1.0203 . 

I.e., the Gompertz Curve is Y = (1.0203) (0.0867) (0*6175) 
Using the 99th percentile, the maximum demand is 





X 

max 


11.3964 + (5-1) 


— 


15.3964 










- 


15. 








TABLE II 








COMPUTATIONS 


IN FINDING A GOMPERTZ 


CURVE 




X 


Y 


Log Y 




I Log Y 




5 


0.1 


- 1.0000 








6 


0.2 


-0.6990 




E 1 log Y = 


-2.0969 


7 


0.4 


-0.3979 








8 


0.6 


-0.2218 1 








9 


0.7 


-0.1549 




E 2 log Y = 


-0.4737 


10 


0.8 


-0.0969 








11 


0.9 


-0.0458 








12 


0.9 


-0.0458 




Z 3 log Y = 


-0.0915 


13 


1.0 


0 









A program for the TI-59 to determine the maximum demand 
when the smallest observed X is unity, and to determine the 
parameters (k, a and b) , is given in Appendix B (Program Bl) . 
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Experimentation with simulated data seems to indicate that 
the Gompertz Curve procedure works well when supply is un- 
limited, but not nearly so well when supply is limited. 

B. USING THE SCHMEISER-DEUTSCH DISTRIBUTION TO ESTIMATE THE 

MAXIMUM DEMAND 

K. Pearson pointed out that a theoretical distribution must 
have at least four free parameters, in order to fit the four 
characteristics (location, dispersion, skewness, and kurtosis) 
adequately [Ref. 6]. Kottas and Lau recommends the four- 
parameter Schmeiser-Deutsch Distribution as most suitable 
[Ref. 6] . The details of the Schmeiser-Deutsch Distribution 
are given in Appendix A. 

The Schmeiser-Deutsch Distribution has a finite upper 
tail. Furthermore, if X is a random variable having the 
Schmeiser-Deutsch Distribution, there exists a closed-form 
expression for the maximum value of X. This expression is 
a+b(l-d) , where a, b, c and d are the four free parameters. 

Using the data in the previous example, we will employ 
the Schmeiser-Deutsch Distribution to estimate the maximum 
demand. We first estimate the parameters using the algorithm 
in Appendix A: 

1) The mode m = (7+8)/2 = 7.5. 

Set a = 7.5. 

2) Set d = (Y(7)+Y(8) )/2 = 0.5. 

3) Select the two desired quantiles which satisfy equation 
(9) , namely, 
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p ^ = 0.2, = 6 

and 

p 2 = 0.9 , x 2 = 11 

Substitute in equation (10), yielding c = 2.9453. 

4) Substitute in equation (11), yielding b = 52.0120. 

Then, the maximum demand is 

= a + b(l-d) C 

max 

= 14.2529 - 14. 

This compares very well with the result obtained by using 

the Gompertz Curve which was X =15.3964. For these data, 
r max 

there is no substantial difference between these two alterna- 
tive methods of estimating X 

3 max 

A program for the TI-59 to determine the maximum demand 
and the four parameters of the Schmeiser-Deutsch Distribution 
is given in Appendix B (Program B2). 

C. ESTIMATING MAXIMUM DEMAND WHEN THERE IS NO DATA ON 
UNSATISFIED DEMAND 

In the previous work it was implicitly assumed that the 
value of demand for each period was obtainable even if demand 
is greater than supply. This rarely happens in the real 
world unless backorders are taken. Speaking of the "newsboy 
problem", the newsboy can collect the data in the range of 
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demand from zero to ordering quantity (Q) only. When the 
newspapers are sold out, he will not know how many more were 
in demand for that day. If he is recording the sales, his 
observed data on that day will be X = Q even if the demand 
is actually X >_ Q. 

The Schmeiser-Deutsch Distribution can be easily used to 
investigate the case when data for demands greater than the 
quantity Q are not available. This is because its parameters 
can be estimated using only the mode and two selected distinct 
percentiles . 

As an example, suppose that the newsboy orders only 10 
newspapers for each day. His sales data for a 10 day period 
are 10, 6, 9, 7, 5, 10, 10, 7, 8, 8. (Instead of demand, 
which was 10, 6, 9, 7, 5, 13, 11, 7, 8, 8.) The sales data 
and their cumulative relative frequency will be: 



TABLE III 



SALES DATA AND THEIR CUMULATIVE RELATIVE FREQUENCY 



X 



Y(X) 



5 



0.1 



6 



0.2 



7 



0.4 



8 



0.6 



9 



0.7 



10 



1.0 
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Since Pr(X = 10) is actually Pr (X >_ 10) , we will ignore 
this data point. 

Following the algorithm in Appendix A for the rest of the 
data points: 

1) Set a = 7.5. 

2) Set d = 0.5. 

3) Select Pj_ = 0.6, = 8, P 2 = 0.7, x 2 = 9 which 

satisfy equation (9) . 

Substitute in equation (10), yielding c = 1.5850. 

4) Substitute in equation (11), yielding b = 19.2279. 

From this, the maximum demand is X = 13.9093 which is 

max 

close to X when the newsboy has unlimited supply, 
max 

With this X , one can determine Q*. Procedures to do 
max 

this are discussed in the next chapter. 
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Ill . ESTIMATION OF OPTIMAL ORDER QUANTITY 



In the previous chapter it was shown how to estimate the 

maximum demand (X ) . In this chapter, the estimate of the 

rricix 

optimal order quantity (Q*) for the case of uncertainty and 
also for the case of risk will be discussed. 

A. OPTIMAL ORDER QUANTITY UNDER UNCERTAINTY 

When data in hand is not sufficient to estimate the 
possible future values of demand or when the distribution 
of demand is not known, the problem of deciding on Q falls 
into the class called decisions under uncertainty. 

It has been shown that the optimal decisions under uncer- 
tainty (Q*) are [Ref. 1]: 



and 



Q* 




X 



max 



for a continuous 
demand , 



* 

Q < 




) (X 



max 



+ 1 ) 



< Q* + 1 , 



for a discrete demand, where c is the unit cost of surplus 
and c q is the unit cost of outage. 

Since one consequence of each decision is another obser- 
vation of demand, one begins under uncertainty and uses the 
uncertainty optimal order quantity until demand data is ade- 
quate to estimate the distribution of demand and change to 
the optimum order quantity under risk. 
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B. OPTIMAL ORDER QUANTITY UNDER RISK 

When one says "under risk", one means that the probability 
distribution of demand is known or, at least, some estimate 
of the probability distribution of demand is available. The 
latter case is possible using the Schmeiser-Deutsch Distri- 
bution suggested in Appendix A. 

The optimal order quantity under risk is the value of Q* 
such that 



F(Q*) 




for a continuous demand 
distribution , 



and 

c 

F(Q*-1) < c "£c" < F(Q * } ' 

s o 

for a discrete demand distribution. 

It should be noted here that all one really needs to 

estimate the probability distribution of demand is a value 

th. 

for the [ c 0 /( c 0 +c s )l quantile. The inverse c.d.f. of the 
Schmeiser-Deutsch Distribution can easily take care of this 
problem. 

As an example, suppose that the newsboy purchased the 

newspapers at the cost of IOC each and he can sell them for 

20C each, with leftovers having a salvage value of 5t each. 

Let p be the ratio [c /(c +c ) ] . Therefore, p = 2/3. 

o o s 

Consider the data from the previous example at the end 
of Chapter II. These data are adequate for a decision under 
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risk, 
c.d.f . 



^ h 

The p quantile can be estimated by the inverse 
of the Schmeiser-Deutsch distribution. 



F -1 (p) = 



a-b (d-p) 



a+b (p-d) 



if p <_ d 
if p > d , 



where the four parameters were found to be a = 7.5, 
b = 19.2279, c = 1.5850 and d = 0.5. 

Therefore, F ^(p) = 8.623, which is the value of Q* for 
the case of continuous demand. 

For the case of discrete demand, the optimal order quan- 
tity under risk is the value of Q* such that F(Q*-1) < p < F (Q* ) , 
and then Q* = 9 in this example. 

In the next chapter we will give conclusions and applica- 
tions from this work. 
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IV. CONCLUSIONS AND APPLICATIONS 



One of the important requirements needed to obtain the 
optimal order quantity in newsboy problems is the nature of 
the demand distribution. This thesis suggested some proce- 
dures which can be used to estimate the demand distribution 
even if data on unsatisfied demand are not available. In 
the case where it is hard to tell what common distribution 
the demand should belong to, the Schmeiser-Deutsch Distribu- 
tion was felt to be useful. 

It is fitting at this point to explain how the procedures 
from this and other research on the newsboy problems would 
be used. 

Suppose we are starting a sequence of single-period 
inventory ordering decisions but have no information on the 
distribution of demand. 

One way of employing these results would involve three 
phases : 

1) For the decisions at time period 1 through some period 

m, the newsboy would probably rely upon a guess at the value 

of the maximum demand (X max ) . Considering the unit cost of 

surplus (c ) and the unit cost of outage (c ) , he has two 
s o 

alternatives : 

(a) In case of c << c , use Q* = X 

so max 

(b) Otherwise, use the uncertainty decision rule. 



)* = 



[c o /(c s 



+ c o )] 



X 



max 
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2) For time periods (m+1) through some later period n, 

he would use Q* = [c /(c +c )] *X , with the value of X 

o o s max max 

obtained from the Gompertz Curve or the Schmeiser-Deutsch 
Distribution, using data from the first m periods. Clearly, 
he may wish to revise this estimate as more data become 
available . 

3) After time period n, he can switch from treating the 
decision as a decision under uncertainty to a decision under 
risk, using the inverse c.d.f. of the Schmeiser-Deutsch Dis- 
tribution to estimate the value of [c / (c Q + c g ) ] quantile, 
which is known to be Q* for continuous demand. 

It should be noted here that because of the data needed, 
m must provide at least three periods where he did not sell 
out, and from the work of Yong-u Sok, n should provide at 
least five periods where he did not sell out [Ref. 4] . 

From the work outlined above we have a reasonable means 
of starting a sequence of newsboy decisions when there is no 
information about the demand distribution and lost sales 
are unknown. It is hoped that the procedures suggested in 
this paper will be useful to those interested in single- 
period inventory problems and their solutions. 
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APPENDIX A 



THE SCHMEISER-DEUTSCH DISTRIBUTION 

In this appendix, the Schmeiser-Deutsch Distribution 
will be briefly reviewed. For a further study [Ref. 6] and 
[Ref. 7] are suggested. 

Due to the lack of convenient methods of generating random 
values from distributions having more arbitrary shapes, the 
Schmeiser-Deutsch Distribution was developed and found suit- 
able for generating random values from many common statistical 
distributions . 

The Schmeiser-Deutsch Distribution can take on a wide 
variety of shapes, which makes it useful for many applica- 
tions. Distributions ranging from Bernoulli trials through 
U-shaped distributions to the uniform distribution, as well 
as heavier tailed distributions, are obtainable. Also, shapes 
ranging from symmetry to maximum skewness, as measured by the 
third moment, are possible. The Schmeiser-Deutsch Distribu- 
tion can have values for the first four moments so as to 
match any of the well-known distributions, and straight- 
forward parameter determination techniques exist. 

A most favorable characteristic of the Schmeiser-Deutsch 
Distributions is the existence of simple closed-form expres- 
sions for both the c.d.f. and the inverse c.d.f., being 
respectively. 
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d- [ (a-x) /b] if (a-bd c ) < x < a 



F (x) 



P 



(7) 



d+ [ (x-a) /b] ^ //c if a <_ x <_ a+b(l-d) C 



and 



a-b(d-p) c if p < d 



x 




( 8 ) 



a+b(p-d) C if p > d 



where a, b, c and d are the four parameters of the distribu- 
tion, and p is a percentile [Ref. 7]. 

"The location and the spread of the distribution are 
determined by a and b, respectively. The shape of the dis- 
tribution, often measured in terms of skewness and tail- 

Q 

weight (or peakness) , is determined by -(d-p) in case of 
p <_ d, or (p-d) c in case of p > d. This factor is a trans- 
formation from the rectangular shape of the uniform distribu- 
tion to the shape of interest. The d = 0.5 yields a symmetric 
distribution. If c > 1, then skew is to the right for 
d < 0.5 and to the left for d>0.5. (If c < 1, the direction 
of skew is reversed.)" [Ref. 7] 

"The distribution possesses two characteristics which may 
adversely affect its usefulness for specific modeling pur- 
poses. The first characteristic is that the value of the 
density function evaluated at the mode, f(x = a), assumes 
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only three values: zero for 0 <_ c < 1, one for c = 1, 

and infinity for c > 1. The second characteristic is the 
truncated tails of the distribution." [Ref. 7] 

"Despite the truncated tails, the distribution's versa- 
tility in assuming a wide variety of shapes makes it a 
reasonable model for a wide range of processes. This is due 
to the ability to assume the essential features of a process, 
which often makes the truncation insignificant." [Ref. 7] 

Parameter Determination [Ref. 7] 

The mode and quantiles are used in estimating parameters 
for the Schmeiser-Deutsch Distribution. The Schmeiser-Deutsch 
Distribution is useful in two common situations. The first 
is when full information is not available, but such quantities 
as the "most likely value" and the "minimum" and "maximum" 
values can be estimated. The second situation is when data 
is in histogram form where F(m) is the percentile of the 
mode m. The two desired quantiles, p^ = F(x^) and p 2 = F(x 2 ), 
must satisfy either of the following two expressions: 



(la-x-J < |x 2 -aj and [d-p-J <|p 2 -d|) 

( I a “ x ]_ I < I x 2 ” a I and ! d -P-L I < |P 2 " d l) 



(9) 



Although computations to determine the four parameters 
can be easily handled by computer, much attention has been 
given in this thesis to develop "manual procedures", meaning 
manual computational capability. The TI-59 handheld calculator 
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(characterized by low price and surprising capability and 
compactness) can be considered as a "manual" tool, and a 
program for the TI-59 to determine the four parameters is 
given in Appendix B. 

An algorithm for estimating the parameters is as follows: 
[Ref. 7] 



1) 


Set 


II 

3 












2) 


Set 


d = F (m) 










3) 


Set 




















a-x. 




d-p 








c 


= In 


-L 


/ In 






(10) 








a-x 2 




a-p 2 






4) 


Set 


b = |a-x 1 |/|d-p 1 j 


l c 




(11) 


5) 


Due 


to the existence 


of the inverse c.d.f. 


in closed- 



form, the value of X corresponding to a percentile p of the 

distribution may be found by inserting p into equation (8) . 

Furthermore, the range of X can be calculated easily 

by substituting p = 0 and p = 1 at Step 5) , or obtained 

c c 

directly from equation (7) since X . = a-bd and X = a+b(l-d) . 

1 ^ mm max 

At this point, one should determine whether or not the ob- 
tained distribution is an adequate model by plotting the 
p.d.f. or c.d.f. [Ref. 7], 
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APPENDIX B 



TI-59 Programs 

This appendix presents two TI-59 programs to determine 
the maximum demand and the parameters. The first program 
is using the Gompertz Curve with the restriction that the 
smallest X is unity. The second program is using the 
Schmeiser-Deutsch Distribution with the two desired quantiles 
satisfying equation (9) . 

Program 

1. Program Description 

The program will estimate the maximum demand using 

the Gompertz Curve, and will also determine the parameters 

b X 

(k, a and b) in Y = ka° . 

2. User Instruction 



Step 


Procedure 


Enter 


Press 


Display 


1 


Enter n 


n 


m 


n 


2 


Enter , . . . , Y^ n 


Y i 

• 


CEZSD 

• 


Y i 

• 






• 

Y 

n 


• 

1 


• 

J^log Y 






Y n+1 

* 


dza 

♦ 


Y n+1 

• 








• 

re7S! 


• 

♦ 

2 2 log Y 
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Step Procedure 



3 Compute X 

c max 

4 Review parameters 
(optional) 

k 

a 

b 



Enter 


Press 


Display 


Y 2n+1 


CS75-] 

• 


y 

2n+l 


• 

• 

Y, 

3n 


• 

re75i 


Z 3 log Y 




cm 


X 

max 



E3 


k 


cm 


a 


El 


b 



The maximum demand shown on the display is for the case 
when the smallest observed X is unity. Otherwise, the maxi- 
mum demand should be adjusted by adding X . -1, where X . 
is the smallest observed X. 
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Program B2 



1 . 
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c and d) 

F(x) 

2 . 



Program Description 
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Distribution. Also 
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